URN TO 
LIBRARY 










TECH LIBRARY KAFB, NM 


□133SLS 


NASA TECHNICAL NOTE 



NASA TN D-7382 


Fssffsssn. 

$ 3.00 


„ 73 - 231 74 

Onclas 

HI/12 0« 19 


R \ Y‘ 

¥§. ‘fe.T’/k 


DISPERSION OF TURBOJET 
ENGINE EXHAUST IN FLIGHT 


by James D. Holdeman 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AftONABfiCS ANB SPACE ADMNttTtATtON • WASHINGTON, B. C • AOGBST W3 




1 Report No 

NASA TN D-7382 


4 Titlp and S-t* »<e 

DISPERSION OF TURBOJET ENGINE EXHAUST IN FLIGHT 


& Rrporr wtt 

August 1973 


6 fVtofmmg O 9* ru ;at i on Code 


) Authotiti 

James D. Holdemun 


0 IV fofrrvnf Of fan* /at .on Name and Ackirm 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44 135 


12 SpooKNinfl A gene* Name and Addttu 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


t Peftofrmog O«gani/ation Report Nj 

E-7479 


10 Work Lkvt No 

501-24 


It Contract or Grant No 


13 Type o< Report and Retrod Covered 
Technical Note 


14 Sponsoring Agan.y Code 



The dispersion of the exhaust of turbojet engines into the atmosphere is estimated by using a 
model developed for the mixing of a round Jet with a parallel flow. The analysis is appropriate 
for determining the spread and dilution of the jet exhaust from the engine exit until it is en- 
trained in the aircraft trailing vortices. Chemical reactions are not expected to be important 
and are not included in the flow model. Calculations of the dispersion of the exhaust plumes of 
three aircraft turbojet engines with and without afterburning at typical flight conditions are 
presented. Calculated average concentrations for the exhaust plume from a single engine jet 
.fighter are shown to be in good agreement with measurements made in the aircraft walce during 
flight. 
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DISPERSION OF TURBOJET ENGINE EXHAUST IN FLIGHT 
by James D. Holdeman 
Lewis Research Center 


SUMMARY 

The dispersion of the exhaust of turbojet engines Into the atmosphere Is modeled as 
a round jet mixing with a parallel flow. The analysis is appropriate for determining the 
spread and dilution of the exhaust plume in the region of the aircraft wake from the en- 
gine exit plane until the exhaust is entrained into the aircraft trailing vortices. Chemical 
reactions are not considered since residence times in this region are too short for photo- 
chemical reactions and engine exhaust static temperatures and pressures at altitude 
flight conditions are below levels required for pji y other significant reactions. An eddy 
viscosity model appropriate for the Jet engine exhaust problem was formulated. It is 
used in an integral momentum method to calculate the dispersion of the exhaust plumes of 
both large and small aircraft turbojet engines with and without afterburning at typical 
flight conditions. 

Measurements of carbon monoxide, oxides of nitrogen, and ozone obtained in the 
wake of an F-104 aircraft flying at Mach 0. 8 at an altitude of 11 kilometers are compared 
with exhaust plume concentrations calculated by using the turbulent diffusion model, es- 
timated concentrations at the engine exhaust plane, and experimentally determined am- 
bient concentrations. The calculated average concentrations are shown to be in good 
agreement with values measured in the wake of the source aircraft. 


INTRODUCTION 


A method for estimating the dispersion of turbojet engine exhaust in flight is pre- 
sented. The method is appropriate for determining the dilution of exhaust products and 
the plume size in the region of the aircraft wake from the engine exit plane until the ex- 
haust is entrained in the aircraft trailing vortices. 

Public concern about the possible adverse Impact of a fleet of high-altitude aircraft 
on the environment has motivated interest in aircraft wakes and stratospheric chemical 
kinetics. The significance of chemical reactions in aircraft wakes is expected to be 







dependent on atmospheric condftlons and nn the concentrations of possible reactants. 

The latter, of course, are determined by the engine exhaust emissions and their disper- 
sion into the atmosphere. 

The jet wake behind an aircraft in flight may be considered in three regions as de- 
scribed in reference 1. These are (1) the jet region, which extends from the engine ex- 
haust plane to the downstream location where the exhaust is entrained in the aircraft 
trailing vortices; (2) the vortex region, which extends from the jet region to vortex 
breakup; and (3) the wake dissipation region, which extends bevond vortex breakup. The 
extent of the jet region is dependent on aircraft configuration arid flight speed. For the 
aircraft and conditions examined i^ reference 1, this region was reported to extend to 
distances of the order of 2 kilometers behind the aircraft. 

Chemical reactions are not considered to be important in the jet region. Residence 
times are too short for any significant photochemical reactions. And static temperatures 
and pressures at the engine exhaust at altitude flight conditions, even with afterburning, 
are low enough to preclude any other significant chemical reactions in the jet. This as- 
sumption is supported by chemical kinetics calculations presented in reference 1 for the 
constant -property Jet core of the GE4 turbojet engine at Mach 2.7 and 19.8 kilometers 
altitude with maximum afterburning. Similar unpublished calculations made at NASA 
Lewis Research Center for various altitudes, flight speeds, and power levels typical of 
the J85 afterburning turbojet engine also support this assumption. 

The model considered in reference 1 for the dilution and spreading of the exhaust in 
the jet region is an empirical formulation from experiments on low-speed, incompress- 
ible jets mixing with an annular coaxial stream at nearly common density and temper- 
ature (ref. 2). Since the turbojet engine exhaust is (at least) sonic and jet temperatures 
are of the order of 5 times the ambient temperature with afterburning, the relations of 
reference 2 are not appropriate. Furthermore, these relations do not reduce to the cor- 
rect asymptotic jet spreading and decay relations for distances from the jet exit where 
the free -stream velocity is large with respect to the difference between the jet centerline 
velocity and the free-stream velocity (ref. 3). 

To overcome some of the shortcomings of previous methods, the present investiga- 
tion was undertaken to develop a model which would be appropriate for the initial condi- 
tions presented by the turbojet engine exhaust and which would yield the known asymptotic 
spreading and decay relations at large distances from the jet exhaust. The analysis is an 
integral momentum method for which the dimensionless velocity difference profile shape 
and the shear stress relation must be specified. The eddy viscosity model used in the 
shear stress relation was formulated by assuming that the eddy viscosity is directly pro- 
portional to the jet mass flow rate plus the absolute value of the entrainment mass flow 
rate and that it is inversely proportional to the jet width. 

Calculations of exhaust plume dispersion for both large and small afterburning tur- 
bojet engines at typical flight conditions are presented. Finally, measurements of oxides 
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of nitrogen (NO x ), carbon monoxide (CO), and ozone (Oj) obtained in the wake of an F - 104 
aircraft flying at Mach 0.8 at an altitude of 11 kilometers are presented and compared 
with concentrations calculated by using the diffusion model, estimated exhaust concentra- 

I 

tions, and experimentally determined ambient concentrations. 

i ANALYSIS 



The exhaust of a turbojet engine is modeled as a round jet mixing with a parallel 
flow. This flow is shown schematically in figure 1. Free mixing of ideal gases is as- 
sumed to occur between an axially symmetric jet (subscript j) and an ambient gas stream 
(subscript •»). The two streams are assumed to be chemically inert. The Prandtl and 
Schmidt numbers throughout the flow are assumed to be equal to 1; thus, the conservation 
equations for momentum, total enthalpy, and species concentration are identical. This 
permits the problem to be solved by considering only the continuity and momentum equa- 
tions. The integral momentum method of Donaldson and Gray (ref. 4) for the mixing of 
an axisymmetric jet with a quiescent ambient gas of a different density has been extended 
to the case where the ambient gas may be moving at a velocity u M . 


Integral Relations 

The integral momentum equation evaluated at the half -value radius is 




PV dr + T 5 r 5 


(1) 


where r g is the radius at which the velocity is one-half the sum of the centerline and 
ambient velocities. A second relation, obtained by evaluating the integral momentum 
equation across the entire width of the flow, requires that the net thrust remains constant 
at all downstream locations. That is, 



ujr dr » PjOjfuj 



All symbols are defined in the appendix. 
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Velocity Profiles 
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To provtde a relation between and r, velocity profiles are assumed. Those 
used In reference 4 are adopted here, except that now the exponential relations apply to 
the velocity difference ratio rather than to a velocity ratio. For the core region (see 
fig. 1), 

for r > r t 


and 




for r at 


(3) 


where r t Is the radius of the constant -property core, r 5 is the half-value radius, and 
X « In 2. For the developed region, 



(4) 


where u c is the velocity on the jet centerline at the streamwise location of interest. 

Because the Prandtl and Schmidt numbers are assumed to be equal to 1, the concentra- / 

tion difference, total enthalpy difference, and velocity difference ratios are identical. 

Thus, equation (4) provides the relation for (c - c w )/(c. - c w ) and (h° - h°)/(h? - h° ) in 
addition to (u^ - u„)/(Uj - u^K 1 

Density 

The density can be expressed in terms of the velocity difference ratios through the 
perfect-gae law and the assumption of identical velocity difference, total enthalpy differ- 
ence, and concentration difference distributions. For economy, let U * (u^ - u m )/ 

(Uj - u w ). The density expression Is then 
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where Cp is the specific heat, m is the molecular weight, h is the enthalpy, and 

u ! » (“j - uj. 


, Eddy Viscosity 

A gradient diffusion model for the shear stress is assumed: thus, 


T 5 “ * 


where « is the eddy viscosity. Numerous eddy viscosity models have been postulated in 
the literature. Those of particular relevance to the present problem are 


Prandtl model (cited in ref. 5): 


c » 0. 025pr 5 |u c - u. 


Ferri model (cited in ref. 5): 


Zelazny model (ref. S): 


< i °. °«r {pu) !(pu) c - <#») J 


°* 036 I*** ' p « u «f r dr 


r 5 (1.0 + 0.«|M c -Mj) 


where if is the Mach number. 







The Prandtl model has been very successful In predicting the mean flow In nearly* 
constant-density cases where u, * 0 and In variable-density cases where u M = 0. 
However, this model predicts decreased shear and hence decreased mixing with decreas- 
ing Uj/u^ for all density ratios. This result is consistent with the experimental data of 
reference 2 for mixing of a jet and an ambient stream at nearly common density, but it is 
contrary to the experimental results of Chriss (ref. 6) for the mixing of low -density (hy- 
drogen) jets in air. 

The Ferri model attempts to account for the density ratio effects and does predict 
the correct mixing trends for both the hydrogen-air and air -air experimental results. 
However, this model predicts segregation of the two streams (no mixing) when (pu)j = 
(pM) x . Unfortunately, the initial conditions posed by turbojet engines to flight are often 
near this condition; hence, the Ferri model is not appropriate. 

The Zelazny model is one of the more recent models and has been shown to be effec- 
tive for predicting the experimental results of most of the jet test cases considered dur- 
ing the Langley Working Conference on Free Turbulent Shear Flows (ref. 7). This model 
Is based on the displacement thickness concept previously used in boundary layer anal- 
ysis. However, tlds model, like the Ferri model, is inappropriate for (pu)j w (pu) K . 

The Zelazny model includes an empirical factor for compressibility and an empirical 
function for variation of e with radial position. This latter factor is not included here 
since only the eddy viscosity at the half -value radius is of interest. Compressibility cor- 
rections have been applied by others also; for example, the eddy viscosity model used 
successfully by Donaldson and Gray (ref. 4) is the Prandtl model with an empirical com- 
pressibility correction applied. 

Because of the inapplicability of previous models to the engine exhaust problem, a 
new eddy viscosity model was developed. The displacement thickness concept (eq. (9)) 
suggests that the eddy viscosity is proportional to the entrainment or detrainment of am- 
bient fluid. Since segregation of the streams is not expected even for cases with zero 
net entrainment, it is postulated in this report that the eddy viscosity should be directly 
proportional to the jet mass flow rate plus Lie absolute value of the entrainment mass 
flow rate and inversely proportional to the half-value radius. In addition, Zelazny's 
compressibility factor (eq. (9)) is adopted. The resultant eddy viscosity relation is 
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r 5 (l.0 + 0.6|M c -Mj) 
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The expression analogous to equation (10), which is appropriate to the core region, is 


0.036 


+ | J <P“ X - P~ u ~ lr dr - { Pj u j ' P» u *> ] 

(r 5 + r'jKl. 0 + O.efMj - Mj) 


The relations given provide a closed set of equations for the solution of the spreading 
and dilution of the Jet from a uniform (slug profile) condition at the jet exit plane to any 
desired downstream location. The solution parallels that given by Donaldson and Gray 
(ref. 4) except that, because of the Inclusion of a moving ambient stream, additional in- 
tegral terms appear and many of the constants are different. 

The integrals in equations (1) and (2) can be evaluated in closed form at any stream - 
wise location. For the constant -property core region, r^ is obtained in terms of r^ 
with equation (2), whence the stream wise distance corresponding to r^ is obtained with 
equation (1). The entire solution for the core region is obtainable in closed form. For 
the developed region, r j is obtained in terms of u c with equation (2). The incremental 
changes in streamwlse distance corresponding to selected incremental changes in center- 
line velocity are obtained with equation (1). The streamwlse distance corresponding to 
any desired u c must then be obtained by numerical integration. 

The method described has been used to calculate the mixing for the jet test cases in 
reference 7. For all cases where the jet dynamic pressure was greater than the free r 
stream dynamic prefers, the agreement between calculated and experimental results 
was very good. 


RESULTS AND DISCUSSION 


Dispersion of Jet Exhaust Plumes 


Calculations ol exhaust plume dispersion have been made for both small and large 
aircraft turbojet engines at typical flight conditions. Exhaust conditions for the J85-GE- 
13 were chosen to represent a small afterburning turbojet engine. The estimated exhaust 
conditions for the GE4 engine which was proposed for the Boeing -2707 SST were used to 
represent a large afterburning turbojet engine. Calculations were also made for the dis- 
persion of the exhaust plume of a J7& engine for comparison with flight data obtained In 
the wake of an F-104 aircraft. Both maximum afterburning and maximum power without 
afterburning were considered for each engine and each flight condition. The following 


7 





■ hn- aKin II isem - i n r fa >— IV 


*d * ,. , mTn mna lm '•fimirmt,, ■■*■ < i i.i M i ^ ■ ittA t i h .l i ! 


usA. 


starting conditions were used: flight Mach number, ambient temperature, Jet exhaust 
pressure ratio, jet total temperature, and primary nozzle diameter. An isentroplc ex- 
pansion was assumed from the sonic primary nozzle to obtain the pressure -matched con- 
dition which was assumed to exist at the engine exit. This expansion yielded the jet- to 
free-stream density and velocity ratios and the exhaust jet diameter which provided the 
required initial conditions for the plume dispersion calculation. These initial conditions 
for each of the engines and flight conditions examined are shown in table I. 

Since the shape of the radial distribution is assumed to be invariant downstream of 
the core region and because the velocity difference and concentration difference ratios 
are assumed to be identical, the exiiaust dispersion is completely defined by relations for 
the centerline concentration difference ratio and the half-value radius as a function of 
streamwise distance from the engine exhaust plane. Because Gaussian profiles are as- 
sumed for the radial distributions, an arbitrary choice for the radial extent of the plume 
region roust be made in order to define an average concentration. The jet edge radius 
r e is defined as the radius to the point where the concentration difference ratio is 0. 001 
times the centerline value. Thus, 


r e »3.157r 5 (12) 

and 

. 0. 14s fe - ' - "\ (13) 

\ c j ‘ c « / \ c ] - <W 

In the far downstream region, where p » and u B » (u^ - uj, the net thrust within 
the radius r e is 0. 999 times the initial thrust at the engine exhaust. 

Dilution in the exhaust plume. - The centerline dilution as a function of distance 
downstream from the engine is shown for the J85, GE4, and J79 engines in figures 2(a), 
3(a), and 4(a), re a ctively. The concentration difference ratio is proportional to x~^'^ 
for centerline dilutions greater than 100:1 (concentration difference ratio <0.01) for all 
engines and conditions examined. This is the asymptotic decay expected for ax'cymmet- 
ric coflowing streams at large streamwise distances. The distance for which (c c - c m )/ 
(Cj - c M ) * 1 is the length of the constant -property jet core region. The jet core lengths 
and the distances to 100:1 centerline and 1000:1 average dilutions are tabulated in table i. 

The downstream distance corresponding to a given dilution varies significantly with 
power level. For the flight conditions and engines examined, the distance required to 
achieve a given dilution with maximum afterburning Is from 1/4 to 1/2 of the distance re- 
quired to achieve the same dilution at maximum power without afterburning. The more 
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rapid mixing with afterburning is caused primarily by the larger shear stress calculated 
for tr.e higher thrust jet which occurs with afterburning. 

For the same flight conditions and power levels, the variation in downstream dis- 
tance to a given dilution for different engines would be expected to be directly related to 
the variation in engine exhaust diameter. However, since nominally equal power levels 
yield different exhaust conditions for different engines, engine operating characteristics 
are also significant. For example for Mach 2 flight with maximum afterburning, the dis- 
tance to 100:1 centerline dilution for the GE4 is 4. 5 times the corresponding distance for 
the J8S, while the engine exhaust diameter for the GE4 is 3. 5 times that for the J85. 

The remaining difference is caused by variations in engine exhaust conditions (i.e. , den- 
sity ratio and velocity ratio). 

Because for centerline dilutions greater than 100:1 the concentration difference ratio 
varies as x -2 ^ 3 and since the average concentration difference ratio is 0. 145 times the 
centerline value at any location, the distance to an average 1000:1 dilution is approx- 
imately 75 percent greater than the distance to a 100:1 centerline dilution, as shown in 
table I. Because the regions of the aircraft wake are often discussed in terms of time 
since aircraft passage, where t * x/u w , the times corresponding to the distances for 
100:1 centerline and 1000:1 average dilution are also givet in table I. For all of the con- 
ditions examined, the time after aircraft passage in which a 1000:1 average dilution is 
achieved is less than 15 seconds, as shown in table I. 

Relation between dilution and spreading . - The half -value radius fr * r 5 , where 
(c - cj/(c c - c m i * 0. 5] is related to the centerline concentration difference ratio (veloc- 
ity difference ratio) by the conservation of momentum (eq. (2)). For large distances 
from the exhaust plane, u w » (u^ - u B ), p » p m , and equation (2) reduces to 


2 (^4 

r 5 _ i M 


Since In this region of the flow (c c - c w )/(cj - c m ) ~ x’*^ 3 , equation (14) gives the ex- 
pected result that the halt-value radius is proportional to x 1 ' 3 . The growth of the half- 
value radius with downstream distance for the engines and conditioas examined is shown 
in figures 2(b), 3(b), and 4(b). 

In the far downstream region, the half-value radius of the plume corresponding to 
any specified centerline dilution can be obtained directly, by using equation (14), from 
the initial diameter of the jet and the Initial mass -velocity ratio. 
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Foe most turbojet engines, engine airflow remains nearly constant from maximum 
power without afterburning through maximum afterburning. Thus, equation (14) provides 
the result that the half-value radius corresponding to a chosen centerline dilution does 
not vary significantly with power level, although the downstream distance to the specified 
dilution ratio decreases as power level increases. 

From the relations between the jet edge radius and the half-value radius (eq. (1211 
and the average and centerline concentration difference ratios (eq. (13)1, the relation be- 
tween r„ and (c„„ - c_)/(c, - in the downstream region becomes 

6 «iV «c j 
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The plume radii corresponding to 1000:1 average dilution and the half-value radii corre- 
sponding to 100:1 centerline dilution for the engines and flight conditions examined are 
given in table I. 

Estimated range of plume concentrations for NO x . CO, and Oj. - Table n compares 
estimated plume and background concentrations for afterburning and nonafterburning con- 
ditions for 1000: 1 dilution in the plume. The engine exhaust concentrations expressed in 
parts per million by volume (ppmv) represent the range of values obtained from tests of 
various engines at simulated altitude flight conditions (refs. 8 to 11). The ambient at- 
mospheric concentration* (ref. 12) and diluted concentrations In the plume are expressed 
in parts per billion by volume (ppbv). Ambient ozene levels n ear the lower end of the 
range shown would be expected in the troposphere, while the higher concentrations would 
be expected in the stratosphere. 


Calculated and Measured Pollution Levels in an Aircraft Jet Wake 

As part of a study being conducted to assess the feasibility of locating and sampling 
engine exhaust products in an aircraft wake, a rendezvous flight between the NASA Ames 
Research Center's CV-990 and the NASA Flight Research Center's F-104 was conducted. 
The flight experiment wa« performed at an altitude of approximately 11 kilometers at 
nominal flight speeds of Mach 0.8. The F-104 as it rendezvoused with the CV-9. 0 is 
shown in figure 5(a). Because of favorable weather conditions, a clearly visible contrail 
was produced by the source aircraft, as shown in figure 5(b). This visible contrail was 
used to position the sensing aircraft (CV-990) in the F-104 wake. The separation 
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between the two aircraft was varied from approximately 2 to 10 kilometers during each 
of two encounters, as shown in figure 6(a). For the first encounter, the F-104 was flown 
without afterburning of the J79 engine. During the second encounter, after turnaround, 
the J79 afterburner was operated. 

The CO, N’O x , and Og data obtained during the rendezvous flight with instruments on 
board the CV-990 are shown in figures 6(b) to (d). The decrease in both CO and NO x 
levels as aircraft separation distance was Increased shows the relative mixing and dis- 
persion of the engine exhaust products with distance. The CO data clearly show the 
higher concentration of this species expected dtring afterburner operation. The max- 
imum measured NO x during the afterburning encounter was approximately equal to that 
measured with no afterburning. Tills result would also be expected, since engine altitude 
emission tests have shown that afterburning has only a small effect on the NO x concentra- 
tion (ref. 10). The ozone data do not show any variations which can be related to the 
F-104 wake. 

The oxides of nitrogen were measure^ with a chemiluminescence monitor with a sen- 
sitivity of 10 ppbv. Carbon monoxide was measured with a fluorescent NDI.R instrument 
with a sensitivity of 200 ppbv. Ozone was measured with both an electrochemical total 
oxidant meter and an ultraviolet absorbtlon monitor. The sensitivities of the ozone in- 
struments were, respectively, 5 and 3 ppbv. From manufacturers* specifications, the 
response time of the CO monitor to 90 percent of a step input is about 30 seconds, while 
the response time of the NO„ monitor to 90 percent of a step input is less than 2 minutes. 
However, because of the turbulence in the F-104 wake, the CV-990 was intermittently 
pitched in and out of the visible contrail, and the sampling Instrumentation was subject to 
a constantly varying input signal. Because of this, the time response of the instruments 
can only be approximated. And since a correction is not required for Interpretation of 
the data obtained in the wake, no correction has been applied to the data shown in 
figure 6. 

The instrumentation used was on board the CV-990 for flight evaluation as part of the 
N*'*' obai Air Sampling Program (GASP). This program is discussed in reference 12. 

Uculated dilution and plume spreading for the J79 engine at Mach 0. 8 with and with- 
out afterburning are shown ii: figure 4. In addition to the centerline dilution and half- 
value radius relations, the average dilution and the corresponding jet edge radius r £ 
are shown. Also the half -wingspan of the F-104 and the CV-990 are indicated in fig- 
ure 4(b) for comparison with the plume radius. The concentration of any nonreacting 
species in the ptun.e at any desired downstream distance can be calculated from the dilu- 
tion relation (fig. 4(a)) and known or estimated concentrations at the engine exhaust and 
in the ambient atmosphere. 

For the separation distances between the CV-990 and F-104 during the flight exper- 
iment (fig. 6(al), average NQ> X concentrations were calculated for engine exhaust concen- 
trations of 20 and 40 ppmv without afterburning and 40 and 80 ppmv with afterburning. 
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The ambient NO x level was assumed to be 1 ppbv. Carbon monoxide average concentra- 
tions were calculated for engine exhaust concentrations of 100 and 300 ppmv without 
afterburning and 1000 and 3000 ppmv with afterburning. An ambient CO level of 50 ppbv 
was assumed. The engine exhaust was assumed to contain no ozone. A mean ambient Oj 
level of 32 ppbv was assumed, based on experimental data obtained preceding and follow- 
ing the rendezvous flight experiment. The average concentrations of CO, NO x , and Og 
calculated for the conditions of the experiment are shown in figures 6(b) to (d). The sep- 
aration distance - flight time relation has been used to present the calculated concentra- 
tions as a function of time for comparison with the measured emission levels in the 
plume. Because the CV-990 was traversing the plume in an unknown pattern, the gas 
sampling probe was exposed to a constantly varying input. Also, as discussed pre- 
viously, the sampling instrumentation does not have instantaneous response. Conse- 
quently, concentrations of the order of the centerline concentrations were not expected 
and not observed. The good agreement between the measured emissions and the cal- 
culated average concentrations should not be interpreted as a calibration of the J79 tail- 
pipe emissions, which were not known; nor should the interpretation be made that the 
sample probe yielded an Integrated emission level as defined by equation (13). Rather, 
the comparison is intended to show that concentrations estimated by assuming that the 
exhaust diffuses as a round jet in a parallel flow with no chemical kinetics are of the 
same order of magnitude as the measured concentration levels. 




SUMMARY OF RESULTS 

The dispersion of the exhaust of turbojet engines into the atmosphere was estimated 
by using a model developed for the mixing of a round jet with a parallel flow. The anal- 
ysis is appropriate for determining the spread and dilution of the engine exhaust from the 
engine exit until it is completely entrained in the aircraft trailing vortex pair. In this 
region, chemical reactions are not expected to be important and were not considered in 
the flow model. 

The analysis is an integral momentum method for which the dimensionless local 
velocity difference profile shapes and the shear stress relation must be specified. The 
eddy viscosity model used in the shear stress relation was formulated by assuming that 
the eddy viscosity is directly proportional to the jet mass flow rate plus the absolute 
value of the entrainment mass flow rate and inversely proportional to the jet width. 

Calculations of the dilution of the exhaust for large and small aircraft turbojet en- 
gines with and without afterburning at typical flight conditions are presented. These re- 
sults indicate that a 1000:1 average dilution may be expected in less than 15 seconds after 
aircraft passage. This dilution corresponds to distances from 0.5 to 7 kilometers 
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downstream of the engine exhaust for the engines and conditions examined. The large 
range depends on both engine size and power level. 

Measurements of carbon monoxide, oxides of nitrogen, and ozone made in the wake 
of an F-104 aircraft flying at Mach 0. 8 at an altitude of 11 kilometers at distances from 
2 to 10 kilometers behind the F-104 are presented. Calculated concentrations for the J79 
turbojet engine plume at conditions of the flight experiment are in good agreement with 
the measured pollutant levels. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 31, 1973, 

501-24. 




APPENDIX - SYMBOLS 

Cp specific heat at constant pressure 

c concentration 

h enthalpy 

h° total enthalpy 

M Mach number 

m molecular weight 

r radius 

r e ' Jet edge radius 
Tf jet core radius 

t Tj jet exit radius 

r 5 half-value radius, where (u - u B )/{u c - u B ) « 0. 5 

r ( pu ) 5 where [pu - (pa) a .y[(pn) e - <pu).J «0.5 

u > streamwise velocity 

U 1 V U - 

x streamwise distance 

€ eddy viscosity 

X In 2 

p gas density 

r shear stress 

Subscripts: 
av average 

c centerline 

j jet exit 

<• free stream 

5 evaluated at r j 
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7 AS LI ;. . OSTIAL CCMSmOrC A HD 100. i CIKTEHUNK DILUTION AND 1000:1 AVERAGE DILUTION RESULTS 


Initial conditions 

R.mIia 

E«la* 

«n* 

n*w mat 

au tutor, 

M 

Altitude, 

km 

Pcwr 

ievrt 

Engine extos sat 
radius, 

r i- 

m 

Density 
ratio, 
p. A, 

Velocity 

ratio. 

Mass-velocity 

ratio, 

| 

Jrt core 
length, 

m 

100 1 Centerline dilution 

1000 1 Average dilution 

Downstream 

distance. 

m 

Time from 
aircraft 

sec 

Half.value 

radius, 

r 9 

m 

Downatre#m 

dl*t*nce. 

tn 

Time from 
aircraft 
(manage. 

see 

Jrt edge 
radius, 

V 

m 

JM 


».s 

w 

e.M 

it 

e.4> 

1.10 

i.i 

MO 

0.4*1 

i.t 

470 

0.1 

1 1 




tM 

M 

.31 

*.• 

.40 

1.41 

3.0 

740 

1.3 

1.4 

1300 

3.3 

1.3 




II.S 

O) 

H 

it 

.41 

lit 

I.I 

390 

47 

I.I 

470 

.3 

I.I 




1 

(kl 

.11 

it 

.*1 

1.44 

3.1 

M0 

I.I 

1.4 

1300 

3.1 

1.4 




1 

(»l 

.11 

1.7 

.11 

1.41 

1.1 

370 

I.I 

1.9 

430 

1.1 

t.l 




* 

fcl 

.11 

3.1 

.11 

.39 

3.7 

410 

3.4 

1.4 

1400 

9.9 

9.3 

GW 

i 


tu 

(*l 

1. 11 

3.9 

0.11 

1. 00 

1.3 

1730 

3.1 

9.9 

3900 

3.9 

39 


i 


1 

<w 

.13 

M 

.99 

1.4S 

IS. « 

3710 

4.7 

i.i 

9400 

9.0 

39 


i 


I 

0) 

.91 

7.4 

.40 

1.90 

9.9 

1*40 

3.1 

9.7 

3100 

3.9 

31 


i 


* 

<w 

.77 

1.1 

.it 

1.31 

U.l 

3930 

4.4 

9.9 

4400 

7 9 

31 

JM 

« 


II.S 

(*l 

e.M 

I.I 

0.10 

1.33 

4.9 

440 

t.l 

3.1 

910 


10.4 




11.0 

w 

.33 

3.0 

.» 

.97 

9.3 

1940 

7.9 

3.0 

3100 

13.3 

11.3 


*MaaIwu« aftartowitac. 
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TABLE n. - PLUME AND AMBIENT CONCENTRATIONS OF OXIDES OF 



f 

X 

( 


I 

1 

i 


% 

I 

• » 





i 

i 


NITROGEN, CARBON MONOXIDE, AND OZONE 


Constituent 

Estimated ( 
atmospheri.c : 
concentration, 
ppbv* 

Estimated plume concentrations 

Without afterburning 

With afterburning 

Engine exhaust, 
b 

ppmv 

1000:1 dilution, 
ppbv 

Engine exhaust, 
ppmv 

1000:1 dilution, 
ppbv 

NO x 

1 to 3 

20 to 100 

20 to 100 

40 to 200 

40 to 200 

CO 

10 to 200 

SO to 300 

60 to 500 

500 to 5000 

500 to 5000 

°3 

10 to 200 

0 

10 to 200 

0 

10 to 200 


* Parts per billion by volume. 
b Part« per million by volume. 


DevriopM riqioo: 
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